O 
O 

(N 

> 
O 

in 



6 

o 

Oh 



Temperature, precipitation and extreme events during the last 

century in Italy 

M. Brunetti 1 , L. Buffoni 2 , F. Mangianti 3 , M. Maugeri 4 , and T. Nanni 1 

1 ISAO-CNR - via Gobetti 101, 1-40129 Bologna - Italy 

2 Osservatorio Astronomico di Brera - via Brera, 28 1-20121 Milan - Italy 

3 Ufiicio Centrale di Ecologia Agraria - via del Caravita, 7A, 1-00186 Roma - Italy 

4 Istituto di Fisica Generate Applicata - via Brera, 28, 1-20121 Milano - Italy 



Abstract. The paper summarises activities within a broad-based research program for the reconstruction of the 
evolution of Italian climate in the twentieth century. The main result of the program is that Italian climate is 
becoming warmer and drier with an increase of both heavy precipitation events and long dry spells. Most of the 
observed signals appear to be due to changes in atmospheric circulation causing an increase in the frequency of 
subtropical anticyclones over the western Mediterranean basin. This hypothesis is also supported by the evolution 
of Italian total cloud amount in the 1951-1996 period. 
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O 1. Introduction 



Around the mid 1990s, the authors set up a broad-based research program with the aim of better understanding the evo- 
, lution of Italian climate in the last 100/150 years. The program was developed both within European (UE IMPROVE 
and ALPCLIM projects) and National projects (National Research Council (CNR) project "Reconstruction of the Past 



o 



CZ) , Climate in the Mediterranean area"). At present it is in progress within the "Progetto Finalizzato CLIMAGRI", a 



project of the Italian "Ministero per le Politiche Agricole e Forestall" . Moreover, in the next two years, further activities 



^ will be performed within the research program " Local climate variability in relation to global climatic change phenom- 
^ (— i ena" funded by the Italian "Ministero per l'lstruzione, l'Universit e la Ricerca" and by Genoa, Milan, Trieste, Turin 

O { . and Udine Universities. The studies so far carried out have improved the availability and the quality of Italian data 
and have produced interesting information on the evolution of temperature, precipitation and some other parameters 
.h in the last 100/150 years. The paper summarises the main results obtained within the research program. 

CZ 2. Data and methods 

In Italy there are several historical series, but only a few of them are digitalized and available for climatological studies, 
while most are still on paper archives. With the exception of some well known series (e.g. Milan, Padua and Rome), 
the construction of the most important database, so far available, of long historical series of Italian meteorological data 
started in the 70s in the framework of a CNR Project. In the second part of the 90s this database was updated and 
improved in the framework of the CNR Special Project " Ricostruzione del clima dell'area mediterranea nel passato". 
The series included in the database are listed in table |l|, and their geographical distribution is shown in figure |l|. All 
series include monthly mean values of daily maximum (Tmax), mean (T) and minimum (Tmin) temperatures and 
monthly total precipitation (P), part of them include daily observations too. 

Monthly series were divided in two groups corresponding to two climatically homogeneous areas - Northern Italy 



(N) and Central-Southern Italy (S) - that are, respectively, the continental and the peninsular zones of Italy (Lc 



Vecchio and Nanni, 1995; Buffoni ct al., 199£). After establishing the new database, the Craddock homogeneity test 
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Table 1. List of the historical series, co-ordinates of the stations and data availability. CNR data base is the result of 
the CNR project "Reconstruction of the past climate in the Mediterranean area". The four last columns of the table 
indicate the goals to be performed within 2003 for monthly and daily series. 







Lat 


Lon 


High 


Monthly CNR database 


Objectives 


1 emperature 


... 
Precipitation 


Monthly data 


Daib 


Temperature 


Precipitation 


Temperature 


1 


Alexandria 


44° 41' 


08° 45' 


98 




1857-1985 


1878-2000 


1875-2000 


1878-2000 


2 


Arezzo 


43° 27' 


11° 52' 


274 


1877-1992 


1876-1992 


1873-2000 


1876-2000 


1879-2000 


3 


Belluno 


46° 07' 


12° 13' 


404 


1873-1993 


1875-1993 


1873-2000 


1875-2000 


1879-2000 


4 


Benevento 


41° 08' 


14° 45' 


135 




1870-1996 




1870-2000 




5 


Bologna 


44° 29' 


11° 20' 


60 


1866-1996 


1813-1988 


1814-2000 


1813-2000 


1814-2000 


6 


Bolzano 


46° 31' 


11° 21' 


272 






1850-2000 


1850-2000 




7 


Bra 


44° 42' 


07° 52' 


308 








1861-2000 




8 


Brixen 


46° 43' 


11° 39' 


569 






1865-2000 


1865-2000 




9 


Cagliari 


39° 12' 


09° 19' 


55 


1879-1996 


1853-1996 


1879-2000 


1853-2000 


1879-2000 


10 


Catania 


37° 30' 


15° 05' 


75 






1901-2000 


1901-2000 


1901-2000 


11 


Cuneo 


44° 24' 


07° 31' 


536 


1877-1993 


1877-1992 


1877-2000 


1877-2000 


1879-2000 


12 


Domodossola 


46° 07' 


08° 17' 


300 








1872-2000 




13 


Ferrara 


44° 49' 


11° 36' 


15 


1879-1996 


1865-1983 


1879-2000 


1865-2000 


1879-2000 


14 


Florence 


43° 47' 


11° 15' 


51 






1889-2000 


1860-2000 


1889-2000 


15 


Foggia 


41° 28' 


15° 32' 


80 






1901-2000 


1901-2000 


1901-2000 


16 


Genoa 


44° 24' 


09° 08' 


53 


1865-1981 


1833-1980 


1833-2000 


1833-2000 


1833-2000 


17 


L'Aquila 


42° 21' 


13° 24' 


735 


1874-1994 


1874-1994 


1874-2000 


1874-2000 


1879-2000 


18 


Livorno 


43° 33' 


10° 18' 


3 


1866-1994 


1857-1994 


1866-2000 


1857-2000 


1870-2000 


19 


Locarno 


46° 10' 


08° 47' 


379 








1901-2000 




20 


Lugano 


46° 00' 


08° 54' 


276 






1864-2000 


1901-2000 




21 


Mantova 


45° 09' 


10° 45' 


20 


1866-1993 


1840-1995 


1828-2000 


1840-2000 


1828-2000 


22 


Messina 


38° 12' 


15° 30' 


54 


1881-1992 


1866-1992 


1881-2000 


1866-2000 


1881-2000 


23 


Milan 


45° 28' 


09° 11' 


122 


1763-1995 


1764-1996 


1763-2000 


1764-2000 


1778-2000 


24 


Marienberg 


46° 44' 


10° 29' 


1323 






1847-2000 


1847-2000 




25 


Naples 


40° 50' 


14° 15' 


149 


1865-1989 


1821-1969 


1865-2000 


1821-2000 


1870-2000 


26 


Padova 


45° 24' 


11° 52' 


14 


1866-1979 


1752-1987 


1866-2000 


1752-2000 


1877-2000 


27 


Palermo 


38° 06' 


13° 21' 


71 


1866-1992 


1806-1992 


1866-2000 


1797-2000 


1876-2000 


28 


Parma 


44° 48' 


10° 18' 


57 


1866-1993 


1833-1987 


1866-2000 


1833-2000 


1878-2000 


29 


Pavia 


45° 10' 


09° 09' 


75 


1866-1980 


1812-1979 


1866-2000 


1812-2000 


1870-2000 


30 


Perugia 


43° 06' 


12° 23' 


520 


1866-1993 


1811-1993 


1866-2000 


1811-2000 


1876-2000 


31 


Pesaro 


43° 52' 


12° 52' 


11 


1871-1991 


1866-1991 


1871-2000 


1866-2000 


1871-2000 


32 


Piacenza 


45° or 


09° 43' 


50 


1872-1995 


1872-1992 


1872-2000 


1872-2000 


1878-2000 


33 


Reggio Calabria 


38° 06' 


15° 39' 


15 


1879-1996 


1867-1996 


1878-2000 


1867-2000 


1878-2000 


34 


Reggio Emilia 


44° 42' 


10° 37' 


62 


1876-1996 


1867-1996 


1876-2000 


1867-2000 


1879-2000 


35 


Riva 


45° 52' 


10° 45' 


73 






1870-2000 


1870-2000 




36 


Roma 


41° 54' 


12° 28' 


56 


1831-1994 


1782-1996 


1831-2000 


1782-2000 


1870-2000 


37 


Rovereto 


45° or 


11° 27' 


206 






1862-2000 


1862-2000 




38 


Rovigo 


45° 03' 


11° 46' 


9 


1879-1993 


1878-1993 


1879-2000 


1878-2000 


1879-2000 


39 


Sassari 


40° 43' 


08° 36' 


224 


1876-1997 


1876-1996 


1876-2000 


1876-2000 


1876-2000 


40 


Siracusa 


37° 03' 


15° 17' 


23 


1878-1988 


1869-1989 


1878-2000 


1869-2000 


1878-2000 


41 


Taranto 


40° 27' 


17° 18' 


22 






1901-2000 


1901-2000 


1901-2000 


42 


Turin 


45° 03' 


07° 40' 


275 


1866-1995 


1803-1996 


1753-2000 


1802-2000 


1753-2000 


43 


Trento 


46° 04' 


11° 07' 


199 






1816-2000 


1816-2000 




44 


Trieste 


45° 39' 


13° 45' 


67 




1841-1996 


1841-2000 


1841-2000 




45 


Vallombrosa 


43° 43' 


11° 00' 


955 




1872-1989 


1872-2000 


1872-2000 


1872-2000 


46 


Venice 


45° 26' 


12° 15' 


1 




1836-1980 


1900-2000 


1836-2000 


1900-2000 



(Craddock 



and Nanni 



1979 ) was applied to the T min, T, Tmax and P series, using regi onal means as reference se ries ( Mauger: 



1998 : Buffoni ct al. , 1999 ). Some series were then homogenized ( Auei , 1992 ; Bohm , 1992 ) both on the 



basis of the test results and the stations' history (metadata). After homogenization, the tempera ture (precipitation) 
s eries were co mpleted over t h e per iod 1865-1996 (1833-1996) by means of a procedure described in Maugeri and Nanni 
( 1998 ) and in Buffoni ct al. ( |1999| ). With the completed data, monthly mean values of the Daily Temperature Range 
(DTR) were calcula t ed from Tmin and Tm ax series. Following the procedure described in Maugeri and Nanni ( |1998|) ; 
Buffoni ct al] ( |1999D ; |Brunctti ct al.| ( |2000a| Jb|) the T, Tmin, Tmax, DTR and P series were then averaged over N and S 
and seasonal and yearly anomalies and their 5-y running means were calculated. Seasonal and yearly N and S average 
anomalies were analysed with the Mann-Kendall non-parametric test, as described in Sneyers (1990), to look for a 
trend. The slopes of the trends were calculated by least square linear fitting. The Mann-Kendall test was also used 
for a progressive analysis of the series Sneyers (1990). The correlation between seasonal and yearly DTR and seasonal 
and yearly precipitation and mean temperature was also performed. Some secular precipitation series were available 
with daily resolution too. For these series the proportion of daily precipitation falling in 5 precipitation class intervals 
was calculated for each year and each season. The resulting series were then analysed for trends, giving particular 
emphasis to heavy and extreme events Brunetti ct al. (2000c). The database resulting from the CNR Special Project 
" Ricostruzione del clima dell'area mediterranea nel passato" should be extended to include more series and further 
metadata. This extension is currently underway in the framework of the " Progetto Finalizzato CLIMAGRI" . Other 
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Fig. 1. Geographic distribution of the stations. Black dots refer to the stations with secular series while white squares 
refer to the stations with series available for the last 50 years. Regions (N and S) and subregions (NW, NEN, NES, 
CE, SO, SI) are indicated too. 



extensions will be performed within the national research program "Local climate variability in relation to global 
climatic change phenomena". Concerning precipitation, another database, covering a shorter period (1951-2000), but 
having a higher spatial resolution, was considered; most of the series were extracted from the Italian Air Force (AM) 
data set, and some others come from the Ufficio Centrale di Ecologia Agraria (UCEA), from the Ufficio Idrografico (UI) 
and from some research projects dealing with the recovery of single series. These series were homogenized, completed, 
validated and grouped firstly into two areas (N and S), then in six more restricted areas by principal component analysis 
Brunetti et al. (2001); the mean regional series and their anomalies were calculated for each area. The distribution 



and the grouping of these stations are shown in fi gure The completion of the series and the regional means were 
performed by already known methods (Karl et al., 1995| , 1998; Brunetti et al., |2001 ) and by a new method discussed 



Brunetti et al 



d2002; ). This new method allows to obtain regional series suitable for the calculation of statistics 
(e.g. drought lengths and frequencies) that can not be obtained using station records completed by methods based 
on random generated data. For the 1951-2000 period, also monthly cloud cover series, extracted from the Italian Air 
Force (AM) data set, were validated, homogenized and analysed for trends. 



3. Results 

3.1. Monthly series trends 

As far as mean temperature is concerned, on a yearly basis there is a positive trend with a 0.99 significance level (si) 
both for N and S; on a seasonal basis, considering a 0.99 si, T has a positive trend in all four seasons in S, while in 
N it has a positive trend in autumn, winter and, using a 0.95 si, in spring. The T trends in the annual temperature 
series, calculated by least squares linear fitting, range from 0.4 C/lOOy for N to 0.7 C/lOOy for S. For the winter season 
the slopes are greater, ranging from 0.7 C/lOOy (N) to 0.9 C/lOOy (S), while for the summer season they are lower 
and in some cases not significant. Even if the data set has been checked for homogeneity and some series have been 
homogenised, trend estimates seem to be critically influenced by data homogeneity. The main problem is that data 
homogenisation performed on a national basis hardly allows inhomogeneities due to changes in national standards to 
be detected, as they concern most of the series in a rather short period. The question is still open, as more metadata 
have to be recovered before giving a conclusive answer. The information available at present suggests that the principal 
inhomogeneities seem to concern the end of the nineteenth century when in many observatories the meteorological 
screens were moved from north facing windows to more open positions and the end of World War II. As discussed in 
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Bohm et al. (2001) the effects of these inhomogeneities seem to cause an underestimation of temperature trends. The 
progressive application of the Mann-Kendall test allows a more detailed analysis of the temperature series trends. A 
complete discussion of this analysis is reported in Maugeri and Nanni (1998); the synthesis of the results is that both 
for N and S the positive temperature trend seems to start around 1920. After 1920 the temperature rises rapidly till 
1950, then it is more or less constant from 1950 to 1985, with only a slight drop in the period 1970-1980. After 1985 
it begins to rise again in all seasons. Concerning DTR, the results of the analysis are discussed in detail in Brunetti 
et al. (pOOOb ). The results of the Mann-Kendall test indicate that DTR has a positive trend (si > 95%) with the only 
exception of winter in N (negative) and of spring and summer in S (not significant). The increase in the annual DTR 
in the period 1865-1996 is weak but significant (0.22 C for N and 0.12 C for S) due to the stronger increase of Tmax 
compared to Tmin. The comparison of these results with the literature shows that the Italian situation is anomalous, 
because generally the DTR is characterised by a negative trend (Karl ct al., 1993| ; Easterling et al, 1997). A detailed 
comparison of the results for the period 1865-1996 is however hampered by the lack of literature data. A more detailed 
analysis of the DTR series can be obtained with the progressive application of the Mann Kendall test. The results 
indicate that in the last decades of the 19th century the DTR trend was generally negative. After the initial decrease, 
in all the seasons (except winter) the DTR trend begins to increase from a date included in the period 1920-1940 for 
N and in the period 1900-1920 for S. Then it continues to increase till around 1970 in N and around 1950 in S. In 
the last decades the trends are generally constant in N, whereas in S they are constant in autumn and decreasing in 
spring and in summer. In winter there is a tendency towards a negative trend in N and a positive one in S. As far as 
precipitation is concerned, the series cover a longer period than for temperature. However, in order to allow a better 
comparison, hereinafter the trend results are presented only for the common period. On a yearly basis a negative trend 
(si 0.99) is evident both for N and S; on a seasonal basis, there is a negative trend in spring, summer and autumn, 
whereas in winter the trend is not significant (S); or it is positive (N). The slopes of the P yearly series, calculated by 
least squares linear fitting, range between -104 mm/100 y for S and -47 mm/100 y for N, giving estimated decreases 
in the period 1866-1995 of 135 and 61 mm. These values correspond, respectively, to 18% and 7% of S and N yearly 
mean values. Both for N and S, spring and autumn have the steepest trends. 



Table 2. Yearly and seasonal temperature, precipitation and daily temperature range trends for northern and southern 
Italy (period 1867-1996) defined by linear regression coefficient (b) and associated error (<7f,). Bold numbers: significance 
level grater than 99%, italic numbers: significance level grater than 95%. 



b±c b 







T 


P 


DTR 






(°C/100y) 


(mm/lOOy) 


(°C/100y) 


Winter 


N 
S 


0.7±0.1 
0.9±0.1 


8+7 


-0.10 ±0.07 


Spring 


N 


0.3±0.1 


-20±7 


0.38 ± 0.01 


S 


0.4±0.1 


-39±5 




Summer 


N 
S 


0.5±0.1 


-7±6 
-17±4 


0.20 ±0.08 


Autumn 


N 
S 


0.5±0.1 
0.8±0.1 


-28±8 
-33±6 


0.29 ± 0.06 


Year 


N 


0.4±0.1 


-47±17 


0.17 ±0.07 


S 


0.7±0.1 


-104±12 


0.09 ±0.06 



In the same way as temperature series, also precipitation series were studied by means of the progressive application 
of the Mann-Kendall test. A complete discussion of this analysis is reported in Buffoni ct al. (1999); the most interesting 
result regards P in S whose high significant negative trend seems to be mainly caused by a strong precipitation decrease 
in the last 50 years. All the trend results are summarized in table 0. 



3.2. Relationships among monthly T, P and DTR 

The comparison of the behaviour of T, P and DTR for the period 1865-1996 (figure 0) has been deeply analysed by 



Brunetti ct al.| (gOOOhj). 

The correlation between yearly and seasonal DTR and P is always negative and highly significant (> 99%) whereas 
the correlation between DTR and T is positive (significance > 95%) in spring and summer for N and in spring, 
summer and autumn for S. The negative P - DTR correlation is more significant in N than in S whereas the positive 
T - DTR correlation is comparable in the two geographical areas. The correlation among T, P and DTR is mainly due 
to high frequent variability, but the same behaviour that is present for the seasonal and yearly data is evident also on 
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Fig. 2. 5-ycar running means of T (thick line), P (histogram) and DTR (thin line). 



longer time scales. Both the yearly and the secular correlated behaviours are probably caused by the same changes in 
atmospheric circulation, with warm and dry conditions (high T and DTR, low P) being related to an increase of the 
frequency of subtropical anticyclones over the western Mediterranean basin. These results are also supported by the 
evolution of Italian total cloud amount in the 1951-1996 period, which is in very good agreement with the evolution 
of the other parameters. By analysing the data of 35 stations we observed that there is a highly significant negative 
trend in yearly and seasonal average cloud amount all over Italy (Maugcri et al, 2001) (table ||). The slope is highest 
in winter. 



Table 3. Results of the application of the Mann-Kendall test and of least square linear fitting to the regional and 
subregional WD and PI series. In order to allow comparison between the different data, the results are expressed 
as ratios (in percentages) between the linear regression coefficient and the mean values in the 1951-1996 period. 
Bold numbers: significance level greater than 95%; non-bold numbers: significance level greater than 90%; when the 
significance level is lower than 90% only the sign of the slope is given. 







NW 


NEN 


NES 


CE 


SO 


IS 


N 


S 




Winter 


-26+17 




-41±15 


-38±13 


-34+11 


-19+9 


-34+15 


-31+10 
























Spring 


+ 


+ 








+ 


+ 


























Summer 


-26±10 


-18±6 




+ 




+ 


-16+7 


+ 




Autumn 










-23±9 










Year 


-14±7 


-12±5 


-18±6 


-12+6 


-21+5 




-14+5 


-14+5 




Winter 


+ 


+ 








+ 


+ 






Spring 


+ 


+ 


+ 


+12±5 


+ 




+ 


+ 


s 




















1 


Summer 


+ 


+ 


+17+7 


+17±7 






+10+4 


+ 


1 


Autumn 


+ 




+ 


+17±8 


+ 


+ 


+ 


+11±5 




Year 


+11±5 


+ 


+9±4 


+10±4 


+ 




+7±3 


+6+3 



3.3. Some results concerning daily precipitation series 



An analysis of 67 sites over 46 years of daily precipitation records for Italy was undertaken by Brunetti et al. (2001) to 
identify any changes in the characteristics of precipitation that may have occurred. The research concerned seasonal 
and yearly precipitation, number of wet days (WDs) and precipitation intensity (PI), and consisted of studying the 
trends both for the station records and for some different area average series. PI was analysed both as amount of 
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precipitation per wet day and by attributing precipitation to 10 class-intervals (Osborn et al. , 2000) removing the 
influence of variations in the number of WDs to yield changes in the underlying shape of the WD amount distribution. 




30 



20 



10 



-10 
-20 




Fig. 3. Results of the application of the Mann-Kendall test and of least square linear fitting to the regional (N and S) 
10 class-interval annual contributions. The results are expressed as percentage variations compared to the mean value 
over the 1951-1996 period. Black bins: significance level greater than 95%; grey bins: significance level greater than 
90%. 



Table 4. Annual and seasonal northern and southern Italy cloud amount trends. The results are expressed in oktas 
in 50 years. 



oktas/50y 



NORTH SOUTH 



Winter 


-1.3±0.3 


-1.1 ±0.2 


Spring 


-0.8±0.2 


-0.8±0.2 


Summer 


-0.9±0.2 


-0.5±0.2 


Autumn 


-0.7±0.3 


-0.9±0.2 


Year 


-0.9±0.1 


-0.8±0.1 



The principal results, summarized in figure |^ and table |], are as follows: 

The number of WDs in the year has a clear and highly significant negative trend all over Italy. It depends mainly 
on the winter, the season which accounts for about 50% of the N decrease, and about 75% of the S one. 
Besides the reduction in the number of WDs, there is a tendency towards an increase in PI. This increase is globally 
less strong and significant than the decrease in the number of WDs, and it is not concentrated in one specific season. 
It is worth noticing that, in winter, the PI trend is positive, but very weak, in N and negative in S. 

. In N, the increase in PI is mainly due to a strong increase in precipitation, falling into the highest class-interval, 
whereas in S, it depends on a larger part of the WD amount distributions. As far as the lower part of the distribution 
is considered, the situation is more similar, and both for N and S, the negative trend is best shown for the sum of 
the three lower classes. 



To search for a trend of dry periods in Italy for the last 50 years, Brunetti et al. (2002) proposed a new methods to 
manage series with missing data without using random generated numbers and analysed the previous data set enlarged 
to 75 stations and updated to 2000. The stations were clustered as already described for the previous database and 
complete daily regional average series were obtained from the incomplete station records, and analysed for droughts. 
Droughts were identified by means of 2 indicators: the longest dry period and the proportion of dry days. The most 
remarkable result is a systematic increase in winter droughts over all of Italy, especially in the North, mainly due to 
the very dry 1987-1993 period (figure |). 
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Fig. 4. Proportion of winter dry days in N and S. In order to highlight the increase in the last 20 years, the averages 
over 1951-1980 and 1981-2000 (dashed lines) are shown too. 

The reliability of the regional series was checked by computing some basic statistics concerning total precipitation, 
rainy days and precipitation intensity and comparing them with the same statistics computed for regional series 
obtained by station records completed with methods based on random number generators. 



4. Conclusions 



Globally, our results support the hypothesis we advanced in some previous papers ( |Brunctti ct al] , |000|,|) and 
confirm that there are strong differences in the response to the recent global warming in Italy and other Central and 
Northern European areas. In fact, besides some common features (temperature increase, tendency toward an increase 
in precipitation intensity), there are other data (cloud cover, daily temperature range, number of wet days) that behave 



conversely. All these differences are probably due to the strengthening in the North Atlantic Oscillation (Hurrell, 1995 



1996) that has accompanied the recent warming and that has caused an increase and a northerly shift in the westerlies, 



with consequent advection of warm and moist air over large areas of Central and Northern Europe (Hurrell, 1995 



Jones et al., 1997), and more frequent anticyclones over its southern part (Rodwell et al., 1999) 
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